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Abstract 

A new totally algebraic formalism based on general, abstract lad- 
der operators has been proposed. This approach heavily grounds in 
the superoperator formalism of Primas. However it is necessary to 
introduce many improvements in his formalism. In this regard, it 
has been introduced a new set of superoperators featured by their 
algebraic structure. Also, two lemmas and one theorem have been 
developed in order to algebraically reformulate the theory on more 
rigorous grounds. Finally, we have been able to build a coherent and 
self-contained formalism independent on any matricial representation, 
removing in this way the degeneracy problem. 
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1 INTRODUCTION 



The fundamental problem in perturbation theory is the solution of the Schrodinger 
equation 

= , (1) 

for the stationary states ^(x, y, z) of a system where the Hamiltonian H is 
split into an unperturbed Hamiltonian H° and a perturbation V. Traditional 
treatments of the theory lean heavily on the expansion of correction to an 
eigenfuction in terms of a complete set of normalized eigenfunctions of H° 
[1-4]. However, the problem can also be formulated in terms of obtaining an 
effective Hamiltonian M = U HU\ with U a unitary operator. The unitary or 
canonical transformation [5] method originated by Van Vleck [6], has been 
adopted by Primas [7], J0rgensen and Pedersen [8], Mukherjee et al. [9] 
and others [10]. The U operator is unitary in the Van Vleck and Primas' 
formalism and produces a Hermitian effective Hamiltonian. 

Murray [11] and Primas [7], have been able to show that any perturbation 
theory can be formulated in the domain of the Lie algebras, in this case 
generated by H° and V. In that concern, the solution of a perturbation 
problem is closely connected with the solution of commutator equations of a 
given type. Further, using the spectral resolution of H°, Primas was able to 
show that the general solution can be written more adequately with the aid 
of the superoperator algebra. 

In the above scenario, our main aim is to recast the superoperator for- 
malism of Primas in an algebraic form using, to that end, the basic theory 
of ladder operators [12] thus our work will be reduced to prove that formally 
it is always possible to build a realization. In Part 2 of this series, we will 
show how particular realizations will lead us to successfully check the present 
approach to of the perturbation theory (AFOPT, Algebraic Formulation of 
the Operator Perturbation Theory) . 

The above AFOPT avoids the matrix representation, since as it is well 
known in the commonly used treatments, the perturbative series and hence 
the expectation values of H, depend crucially on the orthonormal eigenbase 
of H°. 

The outline of the paper is as follow. The treatment begins with the 
definition of the eigenbase {|n )} of H°. Then, the ladder operators defined 
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in this eigenbase have been presented with their main characteristics. At the 
same time in this Sect. 2 the multilinear operators fj+fj": and ff^ff^ have been 
stated. These operators will serve to establish a resolution of any operator 
belonging to the operator space T, whose base has been given by {|n°)}. In 
Sect. 3, two lemmas and one fundamental theorem to of the AFOPT are 
presented. In Sect. 4, the perturbation operator theory is briefly presented. 
This section is followed by a summary and discussions in Sect. 5 Finally, the 
paper ends up with the mnemonic technique in order to write the commutator 
equations. 

2 FORMALISM 

2.1 LADDER OPERATORS 

The full Hamiltonian H is split into a zero-order Hamiltonian H° and a 
perturbation V 

H = H° + \V ,withAe[0,l] (2) 

Orthonormal eigenkets of H° which belong to the zeroth-order eigenspace 
of energy e° n are denoted by \n°) 

#|n°>=e>°> (3) 

As the perturbation is switched on the zero-order eigenkets \n°) evolves 
into orthonormal perturbed eigenkets \n) of energy e n . 

Some time ago, De la Pena and Montemayor [12-16] have shown that 
given the discrete spectral resolution of a linear and Hermitian operator P, 
it is always possible to construct raising and lowering operators associated 
to that operator. Hence, related to H° we have at our disposal the discrete 
eigenbase {|^°)}, thus we may state with all generality 

fj + = J2 c n\n + l)(n\ (4) 

n 

and 

V- = J2<-i\ n - l )( n \ ( 5 ) 
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From the orthonormality condition it is easy to see that fj + and f\- are 
ladder operators 

77+|A:> = c fc |A; + 1> (6) 

f ] _\k)= cU\k-l) (7) 

Now, since f) + and fj- are adjoint to each other, the eigenbase {\n)} is a 
common eigenbase to both operators 57+17- and 17-57+ 

fj + fj-\n) = \c n \ 2 \n) (8) 

rj-fj+\n) = |c„_i| 2 |n) (9) 

The coefficients c n and c* _ 1 are complex number related to the eigenvalues 
of 17+17- and f)-f} + . 

Furthermore, we assume that the eigenvalue spectrum is bounded from 
below and from above [13,17,18] 

o , o , , o 

^0 <~ £\ <~ ' ' ' <~ £jv 

Therefore 

c-i = c M = 

From Eqs. 2.7 and 2.8 it follows that f}+f)_ differs from f)-f) + . In order to 
have only one kind of expressions, we adopt the normal ordering, by which 
the normal product of a set of raising and lowering operators is defined to 
be the product arranged, so that the raising operators are to the left of the 
lowering operators. 

2.2 SUPEROPERATORS 

Now, in order to build the algebraic formulation to of the perturbation the- 
ory, let us introduce the notion of superoperator [7,18,19]. The superoper- 
ator algebra of all linear operators acting on the wavefunction space 7i, is 
a linear vector space, called operator space T. Just as we define mappings 
T : 7i — > 7i called operators, so we can define mappings r : T — > T called 
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superoperators. Both kinds of mappings are linear mappings. Also, linearity, 
the sum and the product by scalar, of superoperators are defined analogously 
to the definitions for the operators. Then it is clear that the superoperator 
space is again a linear space. The foregoing clarification is relevant for forth- 
coming developments of the theory. Actually, let us look for the connection 
between operators and superoperators in the present algebraic approach to 
the perturbation theory. 

So as to do that, let us consider an operator A of the operator space T, 
we will assume that it is possible to write in normal ordering the following 
expansion 



Where now the a mn coefficients will depend on the explicit form of the 



A = E E a mnV+ V- 



,71 



(10) 



m n 



operator A. It is immediate to write: 



n 



(11) 




Then it is possible to show that 



[H°,f)™f)™] = 



(12) 



if m — n, and 





(14) 



from which the results Eq. 2.11 and Eq. 2.12 follow. 
Then it is feasible to define the following operators 



A\\ = E a mm fi™fi™ 



(15) 



m 



and 
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Therefore 



4l = ££<W7?€ ( 16 ) 



A = A\\+A ± (17) 



The operators Ay and A±_ are referred to as the parallel and orthogonal 
components of the operator A relative to H°. They satisfy the next relations: 

[H°,A\\]=6 (18) 

and 

[H°,A ± ]?b (19) 

Since A is any operator belonging to space 7", we have split the operator 
space T into two subspaces 7y and 71 . Where 7jj contains all the operators 
that commute with 77°, and 71 all the operators that do not commute with 
H°. It is necessary to remark that 

Tjj U 71 = T (20) 

and 



7jj n 71 = {0} (21) 

As it has been pointed out, the operator space T is a vector space, there- 
fore Eq. 2.16 may be interpreted as the resolution of operator A into two 
components: one parallel component relative to H° and other orthogonal 
component relative to H°. The above remark contains the key which will 
lead us to prove the theorem about the existence and uniqueness of the in- 
verse of a superoperator T ( see Sect. 3 ). The partitioning that has been 
performed is equivalent to the partitioning in block-diagonal and off-diagonal 
of Primas [7] , and this in turn is the same partitioning as the even and odd 
one of J0rgensen and Pedersen [8]. 
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3 TWO LEMMAS AND ONE THEOREM 



As was distinguished by Murray [20] and by Primas [7], the solution of a 
perturbation problem may be formulated in terms of the solution of the 
commutator equation of the type 

[H°,X} = Y (22) 

where H° is the unperturbed Hamiltonian, Y an operator or function of 
operators and X is an unknown operator that has to be determined. Using 
the spectral resolution of H°, Primas [7] has been able to state the general 
solution for Eq. 3.1 in the language of superoperator, as given by 

X-U(X) = T- 1 (Y) (23) 

In Eq. 3.2 II represents the superoperator that projects from any opera- 
tor, that part which commutes with H°, and T _1 denotes the inverse of the 
superoperator T called derivation superoperator generated by H° [7]. Our 
task will be to reformulate Eq. 3.2 in the abstract ladder operator language. 
If we are able to represent the II, T and superoperators in terms of the 
abstract fj + and fj- ladder operators of the Sect. 2, we will have achieved 
the main goal of the present work. To do that, we would like to state two 
lemmas. Before doing that, we will define II (X) as the parallel projection of 
the X operator. [] 

Definition: For any linear and Hermitian operator X G T the parallel 
projection will be defined by 

U(X) = Y,(n°\X\n )\n )(n°\ (24) 

n 

Lemma 1: Given the abstract ladder operators fj + and r)_ the parallel 
projection superoperator defined over the multilinear operators fj+fj™, m,n = 
0, 1, 2, • • • satisfies the following relation 

m+V-) = SmnVTV- (25) 

Proof: The action of the multilinear operator ff+ffh. on any ket \k) may 
be represented by 

l lL{X), T{X) and r _1 (X) in our notation correspond to (X), k(X) and \{X) in that 
of Primas [7] . 
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f}+ft\k) = X(k;m,n)\k + m-n) 

where X(k; m, n) is a multiplicative factor depending on the powers m 
and n and the quantum number k. By definition 

H(f,™fi n _) = Y,(k\v r +V-\k)\k)(k\ 
k 

and rearranging 

n(f}™f} n _) = ^2\(k;m,n)(k\k + m-n)\k)(k\ 
k 



n(r)™r)-) = T, X (k;m,n)5 mn \k)(k\ 
k 



m+V-)=S mn T l Hk-,rn,n)\k){k\ 
k 



which proves Lemma 1 

The next property derives from the definiton of LT itself : 

U(aA + 0B) = aU(A) + 011(B) (26) 
From Eq. 3.5 and Lemma 1 it is easy to obtain the properties 

n(i) = i,| (27) 

n(i||) = i|| (28) 

n(i ± ) = (29) 
Furthermore, from Eqs.2.16 and 3.4 we may deduce the useful identity 
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A ± = A- n(i) (30) 

Definition : The derivation superoperator T is given by 

T(X) = [H°,X] (31) 

with X eT. 

To study this superoperator, it is necessary to state the following lemma. 

Lemma 2: Given the operator H° and its ladder operators fj + andf}_ the 
derivation superoperator of the multilinear operator 7)™ffi G T satisfies the 
following general form: 

r© n ) = ^"E(C-4)l^fc| (32) 

k 

Proof: By definition of T we get 

T{f,™f, n _)\k) = el +m _ n f,™f£\k) - e° k T%f)»\k) (33) 
Multiplying to the right by the bra (k\ and summing up, it follows 

nvW-)=E(4 + rn-n-4)vW-\k)(k\ 
k 

From which Lemma 2 has been proved. 

The next properties are easily derived from the definition of the Y super- 
operator. 

Since T is a linear superoperator one has 

T(aA + (3B) = aiT (A) + f3T(B) (34) 
Also, it is immediate that 

r(A,|) = 6 (35) 

r(A ± ) + 6 (36) 

and since T is the superoperator which forms the commutator from any 
operator of T with H°, one gets 
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r(AB) = AT(B) + r(i)B (37) 

The superoperator T obtains its name from its derivative properties. 

Some comments must be deserved to the last two lemmas. Firstly, from 
Eq. 3.4 one realizes that the action of U is independent on the physics of the 
system, since the Hamiltonian has not been considered explicitly. Hence the 
superoperator II simply split the entire operator space into two subspaces 
(orthogonal and parallel). Secondly, Eq. 3.6 points out directly, that the 
action of Y has an explicit dependence on H°, due to the presence of the 
transition energy Ae° = £fc +m _„ — which is also an immediate consequence 
of the definition of T itself. 

One very fundamental question to build a coherent and self contained 
algebraic perturbation theory, is to assure the existence of the superoperator 
r _1 in the Primas' theory. Primas has prevented from demostrating this 
relevant theorem because he considers that the inverse superoperator 
has the whole operator space T as its domain [7]. On the contrary, we will 
show that r _1 exists solely in the orthogonal subspace T±G T. Therefore, 
we aim to discover the proper arguments leading to demostrate the existence 
and uniqueness of inverse superoperator. A subject that we will now study 
in somewhat greater detail. 

THEOREM: The inverse superoperator T -1 exists and it is unique, if 
and only if the domain and the range of the linear mapping associated with 
it, can be adequately restricted to the orthogonal subspace 71c T. 

Proof: Since the superoperator T is a linear mapping, it allows us to 
introduce the kernel of a linear mapping [21] and hence the kernel of the 
superoperator T, which we denote by kerT, and that we define as the set of 
all the operators X G T such that T(X) = 0. 

Having in mind that a linear mapping whose kernel is {0}, is injective 
[21,22], we find that r,defined by 

r : T -> T (38) 

with 

r(x) = [H°,x] 

is not an injective mapping. Really, Eqs 3.14 and 3.15 show that kerT = 
7j| ^ {0}. However, it is possible to redefine the domain and the range of the 
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mapping T to the orthogonal subspace, since U(T(X)) = 6. Thus redefining 
the mapping T by : 

r : T ± - T ± (39) 

with 

r(x) = [/7°,x] 

we succeed in getting kerT = {0}. 

Actually, if we assume that an arbitrary orthogonal operator, A e 71, is 
such that A G kerT, then T(A) = 6. But, we know that T(A) ^ 6 if A e 71, 
then the assumption is false. Hence the unique element of the kerT is 0. In 
other words, T is injective. Otherwise, the image and the range of T are the 
same, so T must be surjective. Therefore, the inverse of the T exists and is 
unique. Hence, by fair means we can now write 

r _1 (r(x)) = r(r _1 (x)) = x (40) 

if and only if 

X €71 (41) 
and the Theorem has been proved. 

Lastly the following properties are evident from r -1 , since the linearity of 
r _1 follows from the linearity of T, 

T-\aA + pB) = aT-\A)+ PT-\B) (42) 

Thus the perturbational problem has been reduced to the finding of an 
explicit expression for In Part 2 of this series, we will study particular 
forms for (also for V and n), depending on the algebra of ladder operators 
associated to the physical problem to be tackled. 



4 PERTURBATION METHOD 

As aforementioned the complete Hamiltonian H has been split into an un- 
perturbed Hamiltonian H° and a perturbation operator V scaling with the 
real parameter AG [0, 1] 
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H = H° + XV (43) 

Besides, the comments that have been made at the begining of Sect. 2 
(cf. Eqs. 2.1 and 2.2) also special mention deserves the fact that in general 

[H°,V}^0 (44) 

which implies that we cannot find a common eigenbase for H° and V. 
But we can think of a certain unitary transformation, that will change this 
situation. 

The idea of choosing a unitary transformation corresponds to the need 
of leaving invariant the spectrum of eigenvalues of the energy. The unitary 
transformation only modifies the eigenvectors. 

Let U be a unitary transformation defined as 

UHfp = U(H° + \V)fr (45) 

We can now introduce two new operators M and W, through the defini- 
tions 

M = UHU^ (46) 

and 

W — M — H° (47) 
The relation 4.4 allows to write 

M = H° + W (48) 

From Eq. 4.4 it is immediate to see that M has the same spectrum of 
eigenvalues as the Hamiltonian H. 

We will now suppose that U satisfies the following condition 

[H°,W}=6 (49) 

That means that H° and M will have common eigenvectors as follows 
from Eq. 4.6. Therefore, if Eq. 4.7 holds, we may write 

(n°\M\n°) = (n°\H°\n } + (n°\W\n°) (50) 
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e n = e° n + (n°\W\n°) (51) 
Since (n°|M|n°) = e n and M = UHU^ we may write 

UHU ] \n°) = e n \n°) (52) 

Therefore, after multiplying to the left by W and having in mind that U 
is a unitary transformation 

HUi\n°) = e n U ] \n°) (53) 
where W\n°) is the new eigenket of H. 

Briefly, imposing the condition given by Eq. 4.7 we have the following 
scheme: 



In) 



fr\n°) 



(54) 



e n = e° n +(n°\W\n° 



That is to say, resolving the eigenvalue problem for the Hamiltonian H 
implies to find the transformation W that makes possible the Eq. 4.7 which 
in turns, will allow us to write the explicit form of W. 

Let us suppose now that the unitary transformation may be written as 
the exponential of a certain antihermitian operator, G = — G\ henceforth 
referred to as the generator of the transformation. Then we immediately get, 
the relation 

W = M - H° 
W = UHU^ - H° 



W = exp(G)ifexp(-G) - H° 
Using the expansion of Baker-Camppell-Hausdorff 

w=(h + ±[g,A\ + ±[g, [g,h}} + .. 



we get 
H° 



(55) 



(56) 



13 



From Eq. 4.1 we arrive at 

W = XV + H° + AV] + [G, H° + At>]] + • • • (57) 
Let us now assume that 

W = XW 1 + A 2 #2 + • • • (58) 



and 



G = XGi + A 2 G 2 + • • • (59) 



Insertion of Eq. 4.15 and 4.16. in Eq. 4.13, furthermore, developing, 
rearranging and comparing equal powers in A, lead us in a straighforward 
way to 

[H°,G 1 ] = V-W 1 (60) 
[H°, G 2 ] = ±[Gi, V] + I[G ls [G l5 if°]] - W 2 (61) 



[^Gs] = ^[G 2 ,y] + ^[G 1 ,[G 2 ,^]] + i[G 2 ,[G 1 ,^ ]]+ (62) 
i[G l5 [G 1; n] + i[Gx, [G l5 [Gx, iT]]] - W 3 

... an so on. 

It is apparent that the set of last Eqs. 4.18-4.20 is a system of cou- 
pled commutator equations for the G n operators. This set obeys the general 
structure 

[H°,G n ] = A n -W n (63) 

where H° and A\ = V, constitute the data of the problem and the G n 
are the unknown operators to be determined. The A n operators, with 1, 
are specified in terms of H° and A m with m < n. 

It is necessary to determine the W operator, provided that [H°, W\ = 
or equivalently to that of II(W') = W. However, these conditions are fulfilled 
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if, in turn each one of W n results to be a parallel component operator relative 
to H° . On this basis it may be concluded that 

IT(Wg = W n (64) 
Now the operation with II on Eq. 4.21 leads to 

U([H°,G n ]) = U(A n )-U(W n ) (65) 
Having in mind the identity 

U([H°,G n ])=0 (66) 

we get 

n(A0 = u(w n ) (67) 

Thus from Eq. 4.22 we write 

W n = U(A n ) (68) 
Otherwise, from the definition 2 we have that 

T(G n ) = [H°,G n ] (69) 

provided that G n G T±, for every n. However, this condition is equivalent 
to say that 

n(G n ) = 6 (70) 

Therefore from Eq. 4.22 we obtain 

r(G n ) = A n - W n (71) 

or 

r(G n ) = i n - n(i n ) (72) 

But the hand right side of the above equation is an operator that belongs 
to Tl, therefore V is well-defined. Thus, it may be deduced that T" 1 exists, 
in brief 
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r- 1 (r(G n )) = r- 1 (^ n -n(A n )) 



or 



G n = r-\A n - n(A0) 

To sum up, given a problem of the type 



H = H° + V 



we will have that 



and 



e n = e° n + (n°\W\n°) 



Where 



and 



W = XW 1 + X 2 W 2 + 



w n = n(i n ) 



G = Ad + A 2 G 2 + 



G n = r- 1 ^ - n(A0) 

The explicit forms of any A n are : 

A 1 = V 



A 2 = ^[G 1 ,V] + ^[G 1 ,[G ll H ]] 
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A = ^Vj + ^G^H^ + ^iG^H }} (84) 
+^[G 1 ,[G 1 ,V }} + ^[G 1 ,[G 1 ,[G 1 ,H°}}} 

... an so on. 

In order to know all the terms of the series, we have developed a mnemonic 
method (Cf. appendix): 

A = (l) 
A 2 = (i|i)©(i,i|o) 

A 3 = (2|1) © (1, 2|0) © (2, 1|0) © (1, 1|1) © (1, 1, 1|0) 

5 SUMMARY AND DISCUSSIONS 

It has been shown that from the spectral resolution of H°, the abstract 
ladder operators fj + and fj- may be defined. In turn, these operators serve 
to build multilinear operators in normal ordering f)™fft. Taking advantage 
of the properties of the ff^ffl in relation to H°, we have been able to split 
the entire space T into two subspaces 7jj and T± accordingly to any operator 
that commutes or not with H°. The above splitting of T has allowed us to 
demonstrate the existence and uniqueness of T _1 under the condition that 
the domain and the range of T must be the orthogonal subspace T± C T. 
Primas [7] was prevented from demostrating this relevant theorem, because 
he had considered that the superoperator has the whole operator space 
T as its domain. 

As may be seen from Sect. 4, the entire algebraic formulation of the oper- 
ator perturbation method lean heavily on the well-defined 11(77+7/™), T{fj r ^fj r !i) 
and r _1 (r/+r}™) operators. 

As was remarked at the begining, the present approach has been built 
independently on whatever matricial representation. Therefore, the Hamil- 
tonian H° may have any degeneracy, however this situation is immaterial in 
that concern the purely algebraic relations between the operators involved. 
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In Part 2 of this series, the method is seccesfully applied to two quantum 
mechanical systems: "The Stark Effect in the Harmonic Oscillator" and "The 
Generalized Zeeman Effect". 

6 APPENDIX 

In order to write out efficiently the explicit form of the commutator equations 
determining the A n operators, we have developed a mnemonic method. 
RULE 1: A bracket of two sides is drawn 

(•••!•••) 

RULE 2: In the right side we must put 1 or 0. 

RULE 3: In left side of the bracket we must put integers, in such way that 
its sum must be n, i.e. the order of the iteration, consequently the subindex 
of A n superoperator. 

RULE 4: We return to rule 1 until exhausting the possibilities of gener- 
ating further diagrams. 

RULE 5: In order to write an explicit commutator form for each operator 
A n , we must consider 

Left Side 

Besides, we have to remember that each expression is divided by the 
factorial of the number of integers in left side. 
As an example we calculate A 2 and A 3 : 

A 2 = (i|i)e(i,i|o) 

A 2 = ±[G 1 M + ±[G 1 ,[G 1 ,H ]] 
A 3 = (2|1) © (1, 2|0) © (2, 1|0) © (1, 1|1) © (1, 1, 1|0) 



Right Side 

-> H° 

1 -> V 
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As = ±IG 2 ,V] + ±\G U [G 2 ,H°]] + ±[G 2 , [G 1 ,H°]]+ 
l i [G 1 ,[G 1 ,V°]] + ^[G^i^H ]]] 

In what follows we display some diagrams: 

M = (i) 

i 2 = (l|l)©(l,l|0) 
A 3 = (2|1) © (1, 2|0) © (2, 1|0) © (1, 1|1) © (1, 1, 1|0) 
i 4 = (3|i) © (1, 3|0) © (3, 1|0) © (2,2|0) © (1,2|1) © (2, 1|1)© 
(1, 1, 2|0) © (1, 2, 1|0) © (2, 1, 1|0) ©(1,1, 1|1) © (1, 1, 1, 1|0) 
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